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ABSTRACT 



Phenomenological models for evolution of dusty galaxies and E/SO galaxies, 
respectively, are developed to address two major questions concenring galaxy pop- 
ulations in deep infrared (IR) surveys: (1) Do normal late-type galaxies or star- 
burst galaxies (including galaxies with obscured AGNs) dominate among sources 
in deep IR surveys? (2) How much do E/SO galaxies contribute to the counts in 
deep mid-infrared (MIR: 3 - 20//m) surveys? Among three new models for evo- 
lution of dusty galaxies, it is assumed in Model SI that starburst galaxies are the 
dominant population, and in Model S2 that normal galaxies dominate. Model S3 
is an intermediate model. Comparing the model predictions with a wide range 
of observational data collected from the literature, we find that none of these 
models can be ruled out, given the uncertainties of the data. We show that the 
most direct method to distinguish these models is to compare the predicted color 
distributions of IR galaxies with observations, which will soon be available from 
the SWIRE survey. The models for E/SO galaxies follow a simple passive evolu- 
tion approach. Among the three E/SO models (El, E2 and E3) investigated in 
this paper. Model E2 which is specified by a peak formation redshift Zpeak = 2, 
and an e-folding formation time scale to — 2 Gyr, fits the data best. This sug- 
gests a synchronization between the evolution of E/SO galaxies and of starburst 
galaxies, in the sense that the peak of the formation function of E/SOs {zpeak = 2) 
is close to the peak of the evolution functions of starburst galaxies {zpeak — 1-4). 
We find that E/SOs contribute about 10 - 30% of the counts in the MIR bands 
of < lOfim, and up to 30 - 50% of the optical/NIR counts in the bright end. 
Their contributions to counts in the UV (2000A) and in the longer wavelength 
IR (> 12//m) bands are negligible. Taking into account this contribution, new 
predictions for counts and confusion hmits in the SIRTF bands are presented. 
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Subject headings: galaxies: evolution - starburst - elliptical and lenticular - 
Seyfert - luminosity function; infrared: galaxies 

1. Introduction 

Many new windows have recently been opened in various wavebands for observations 
of high redshift galaxies. Particularly, deep ISO surveys (e.g. Rowan-Robinson et al. 1997; 
Pranceschini et al. 1997; Kawara et al. 1998; Elbaz et al. 1999; Aussel et al. 1999; Flores 
et al. 1999b; Puget et al. 1999; Dole et al. 2001; Clements et al. 1999; Serjeant et al. 2000) 
in the infrared and SCUBA surveys (Hughes et al. 1998; Barger et al. 1998; Blain et al. 
1999) in the submm wavebands have shed lights on the dark side of galaxy formation and 
evolution. This stimulated a surging wave of empirical models (Xu et al. 1998; Blain et 
al. 1999; Roche & Bales 1999; Dole et al. 2001; Xu 2000; Rowan-Robinson 2001; Xu et al. 
2001; Pranceschini et al. 2001; Pearson 2001; Malkan & Stecker 2001; Chary & Elbaz 2001) 
to interpret the new ISO and SCUBA sources. The main aim of these empirical models is 
to better constrain some important characteristics of the galaxy evolution, including the 
peak epoch of the star formation rate in the universe and the roles played by different star 
formation modes (e.g. the quiescent mode against the interaction stimulated mode), for 
which the earlier optical deep surveys such as the Hubble Deep Pield Survey (WiUiams et 
al. 1996) may have provided biased information (Madau et al. 1996) by neglecting the dust 
extinction (Rowan- Robinson et al. 1997; Lonsdale 1999). This will certainly have impact 
on the theoretical simulations of galaxy formation and evolution, which (e.g. Somerville et 
al. 2001) have just started facing the fact that most of the early star formation may be 
hidden behind a thick veil of dust, making the incorporation of effects of dust extinction 
and emission in the framework essential (Granato et al. 2000). 

Some of the most recent empirical models (Rowan- Robinson 2001; Xu et al. 2001; 
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Pranceschini et al. 2001; Pearson 2001; Chary & Elbaz 2001) have the feature that they can 
be constrained by counts and the cosmic background radiation in various IR and submm 
bands simultaneously. This is achieved by using SED templates to link sources in different 
bands. The SED library of Xu et al. (2001) contains reahstic SEDs of a complete sample of 
837 IRAS 25//m selected galaxies, enabling prediction of counts as well as color distributions, 
which provide additional constraints to the model. Templates in Rowan- Robinson (2001), 
Xu et al. (2001) and Chary & Elbaz (2001) extend from the IR-submm to optical and 
UV wavebands, hence these models can predict also the contributions from IR sources to 
optical and UV counts, relating the evolution of IR sources to that of galaxies seen in earher 
optical and UV surveys. All of these works found strong evolution among IR sources in the 
redshift range of < z < 1.5. Assuming that the narrow sub-mJy bump on the Euclidian 
normalized differential counts of ISOCAM 15//m surveys (Elbaz et al. 1999) is due to the 
K-corrections caused by the strong unidentified broad band emission features (UIB) at 6 - 
8//m, which are shifted into the ISOCAM 15/xm band filter (LW3) when ^ ~ 1, Xu (2000) 
concluded that typical IR galaxies at z ~ 1 have ^15^^ ~ 10^^ Lq, namely about 20 times 
more luminous than their local counterparts, while their comoving density is about the 
same as their local counterparts. As argued by Xu et al. (2001), the location of the 15/xm 
bump {fi5fj,m ~ OAmJy) provides a strong constraint on the luminosities of IR sources at 
2; ~ 1. This has been confirmed by the detailed study of ISOCAM sources in the HDF-north 
by Elbaz ct al. (2002), who obtained redshifts for nearly all of these sources (40 out of total 
41). The strong evolution has also been supported by surveys in other ISOCAM bands. 
For example, Clements et al. (2001) found that the luminosity functions of their ISO 12//m 
sources are consistent with a pure luminosity evolution of rate ~ {1 + z)^'^, as derived by 
Xu (2000) from the ISOCAM 15/im counts. 

On the other hand, different authors identified different populations of galaxies as 
the major carriers of the evolution of IR sources. It was found in early IRAS studies 
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(Pranceschini et al. 1988; Rowan-Robinson & Crawford 1989) that IR sources can be divided 
into 3 different populations: normal late- type galaxies ('cirrus galaxies'), starburst galaxies, 
and galaxies with AGNs. Assuming that all IR sources undergo pure luminosity evolution 
with the same evolution rate, Rowan-Robinson (2001) found that 'cirrus galaxies' dominate 
the ISOCAM 15//m counts, the SCUBA 850//m counts, and the cosmic IR background 
radiation. Xu et al. (2001), motivated by results of optical identifications of ISO galaxies 
by Flores et al. (1999b) and Aussel et al. (1999) which show a larger percentage of these 
sources are in interacting systems, assumed that most of the evolution of the IR sources 
is due to starburst galaxies (defined by warm FIR colors: fmimj fimum > 0.5). Obscured 
AGNs with relatively low /25//m//60/im ratios (due to the very high extinction affecting 
even the MIR fluxes) may be misclassified as starbursts in Xu et al. (2001), who classified 
galaxies with AGNs using the criterion f25nm/ feoum > 0.2. Pranceschini et al. (2001) also 
assumed that the strong evolution is confined to the starburst population. Their starburst 
galaxies include all 'active' galaxies not classified as Sey 1, i.e. including the Sey 2s and 
the LINERs, in Rush et al. (1993). Pearson (2001) found that a separate population 
of ultra- luminous galaxies (ULIRGs) with L~ IO^^Lq, confined in a very narrow range 
of redshift centered at 2; ~ 1, is mostly responsible for the strong evolution seen in the 
ISOCAM 15//m counts. Since all these models can fit the IR-submm counts and the CIB 
within uncertainty limits, there is apparently a degeneracy concerning the population of 
galaxies that carry most of the evolution of IR sources. The degeneracy still remains even 
when the new, detailed information of the 15 micron universe (source counts, redshift 
distributions, luminosity functions at different redshifts, etc.; Elbaz et al. 2002; Pranceschini 
et al. 2001) is considered (see Section 2). 

One of the central issues in the current agenda of hierarchical galaxy formation 
simulation studies concerns the roles played by two different star formation modes (eg. 
Somerville et al. 2001; Kauffmann et al. 2001), namely the quiescent star formation mode 
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and the interaction-induced star formation mode. There are apparent hnks between the 
normal quiescent galaxies and the quiescent star formation (as in the Milky Way and 
M31) and between starburst galaxies and the interaction-induced star formation (as in the 
Antennae Galaxies and in M82) . Therefore, the answer to the question whether the normal 
galaxies or starburst galaxies dominate the faint IR counts will have important impact on 
galaxy evolution theories. 

Another deficiency of the current models for IR sources is the neglection of early type 
galaxies. Although these galaxies have little dust emission, therefore do not contribute 
significantly to IR counts at wavelengths longer than ~ lO/im, they are found as an 
important population in ISO d.TfjLm counts (Flores et al. 1999b; Aussel et al. 1999; Serjeant 
et al. 2000). Future mid-IR surveys such as those planned for SIRTF IRAC cameras (4 
bands centered at 3.6/im, 4.5/im, 5.8/xm and S/^m, respectively) will certainly detect many 
of these galaxies. 

In this paper, we shall address the above two issues, namely (1) developing new models 
to investigate how to break the degeneracy concerning the major evolution population 
in IR counts (using SIRTF data in particular), and (2) modeling the evolution of E/SO 
galaxies and predicting their contributions to the MIR counts. The plan of the paper is 
as follows: in Section 2 a set of new models for evolution of dusty galaxies are presented, 
while in Section 3 models for evolution of E/SO galaxies are developed and compared to 
the observations; predictions for the UV, optical and NIR counts by a set of composite 
models, each consisting of a dusty galaxy evolution model and an E/SO evolution model, are 
presented and compared with data in Section 4. The predictions for counts and confusion 
hmits in SIRTF bands by the same models are presented in Section 5. In Section 6 we 
show how to use color distributions of ISO sources and future SIRTF sources to answer 
the question whether starburst galaxies or normal late-type galaxies are dominant in the 
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IR sources. Section 7 is devoted to discussion. A summary is given in Section 8.^ The 
A-Cosmology (Ha = 0.7, = 0.3, Ho=75 km sec"-*^ Mpc~^) is assumed throughout the 
paper. 

2. New Models for Evolution of Dusty Galeixies 

2.1. Previous Models 

In all the multi-band models for the evolution of dusty galaxies in the literature, 
the different populations of galaxies are separated by their characteristic spectral energy 
distributions (SEDs): the normal galaxies have relatively low fmfim/ fioofj.m and /25/im//i2/im 
ratios and very prominent UIB features, the starburst galaxies have relatively high 
hoixm/ fioof^m and /25Mm//i2,im ratios and less prominent UIB features over a steep rising 
MIR continuum, and galaxies with AGNs have relatively high f25nm/ feoixm ratios and very 
week UIB (contributed by ISM dust not associated with the AGN). Models by Xu et al. 
(2001) have also considered the luminosity dependence of the SED in each population, while 
some other models (Rowan- Robinson 20001; Pranceschini et al. 2001; Pearson 2001; Chary 
& Elbaz 2001) separate high luminosity starbursts (Arp220 type) from moderate luminosity 
starbursts (M82 type). Because of their different SEDs, these different populations of 
galaxies give different relative contributions to counts at different wavelengths. However, 
it appears that the counts are rather loose constraints to these relative contributions. At 
least 2 factors are responsible for this degeneracy: (1) there are many parameters in the 
evolution functions (both the luminosity evolution and the density evolution) when different 
populations are allowed to evolve differently, and (2) there are still significant uncertainties 

^More detailed results on predicitons of models presented in this paper are available on 
request to C.K.X.. 
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in the current ISO counts (e.g. counts of the same ELAIS 15//m sources reported by 
Serjeant et al. 2000 and by Gruppioni et al. 2002 differ with each other by as much as a 
factor of 3), due both to cahbration errors and the field-to-field variation. In principle, the 
degeneracy can be broken by comparing model predictions on redshift distributions of deep 
FIR (A ^ 60//m) selected samples with observations: Models assuming warmer IR SEDs 
(starbursts dominant) usually predict larger mean redshifts than models assuming cooler 
IR SEDs (normal-galaxies dominant) do. However, similar to the SCUBA galaxies, the 
large beams of FIR detectors (such as the ISOPHOT cameras) make it rather difficult to 
pin-down the optical counterparts (usually quite faint) of faint FIR sources. 

A more direct method to constrain the relative contributions of different populations 
to the IR counts is to compare the predicted and observed color distributions of IR sources. 
Among the multi-band models for evolution of dusty galaxies in the literature, those of Xu 
et al. (2001) have the most sophisticated algorithm in deahng with the SEDs, and are the 
only ones that can predict both the mean colors and their dispersions. Therefore, we choose 
to build our new models using the same algorithm as Xu et al. (2001). 



In the new models, the simulation code of Xu et al. (2001) is modified in the following 
aspects: 

(1) The luminosity and density evolution functions for the star burst galaxies and normal 
late-tpye galaxies, respectively, have the following new form: 



2.2. 



New Models 




(1 + zf^ X ((1 + (1 + z)/{l + Zi))/{1 + 1/(1 + Zi))) 



X 



((1 + 1/(1 + z2))/{i + (1 + z)/{i + Z2)y^^'"^ {z<7y, 



(1) 



-9- 



and 



Gi{z) 



(1 + zy* X ((1 + (1 + z)/{i + zi))/{i + 1/(1 + zi))y^-p* X 



((1 + 1/(1 + Z2))/{1 + (1 + z)/{l + z^)Y^^^ {z < 7). 



(2) 



These smoothly-joined 3-piece power laws, in contrast with the sharply- joined 2-piece 
power laws used in Xu et al. (2001), will allow softening in the low redshift end 
{z < 0.5) of the evolution functions to improve the fit to the redshift distribution of 
IRAS 60//m sources. 

(2) The 25/xm local luminosity functions (LLF) of 3 populations used in Xu et al. (2001) 

did not themselves take into account of the evolutionary effects. In the new models, 
new 25//m LLPs corrected for these effects are used (see Appendix A for the details 
of the new LLPs). It is found that these new LLPs are only marginally different from 
the old ones. 

(3) The UV portion (A < 3000A) of the SEDs in the SED lib, which is an important part 

of the code, is better constrained in this work (Appendix B). 

Given the large number of parameters, it is out the scope of this paper to explore the 
entire parameter space for models allowed by available data. Instead, we concentrate on 
three models which predict very different relative contributions by starburst galaxies and 
normal late-type galaxies to the counts in the IR surveys. Por each of them, parameters 
have been fine-tuned so that the model can fit all available data as well as possible. These 
models for dusty galaxies are presented in Table 1 and Pig.l. 

Model SI is similar to the models in Xu et al. (2001) and in Pranceschini et al. (2001) 
in the sense that the starburst galaxies are assumed to dominate the evolution of the IR 
sources. In Model S2, similar to Rowan-Robinson (2001), all three populations of dusty 
galaxies are assumed to have pure luminosity evolution with the same evolution rate. This 
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is actually the same evolutionary scenario adopted by Xu (2000), who assumed that the 
entire body of IR sources evolve as a single population. Note Rowan-Robinson (2001) 
assumes that there is no SED evolution, while in Model S2 we assume SEDs evolving with 
luminosity. In Model S3, it is assumed that the normal galaxies and star burst galaxies 
give about equal contributions to the ISOCAM 15//m counts, as hinted at by the optical 
identifications of ISO sources (Flores et al. 1999b; Aussel et al. 1999). For Model SI and 
Model S3, the luminosity evolution function of galaxies with AGNs is adopted from the 
optical QSO luminosity evolution function of Boyle et al. (2000) 

Fagn{z) = lO'-''^-'-''^' {z<7). (3) 

This is shghtly different from the power law function used in Xu et al. (2001). 

Tablel. New models for dusty galgixies 

Normals Starbursts AGNs 

model zi zi Ui Vi Wi U2 V2 W2 P2 ?2 r2 Us Vs W3 

51 0.5 0.85 1.5 1.5 2 2 11.5 1.2 1 9.8 1.2 1.36^-0.27^2 

52 0.5 1.2 3.2 5.7 2.5 3.2 5.7 2.5 3.2 5.7 2.5 

53 0.5 1 2 8 2 2 10.2 1.5 2 3 1.5 1.36^-0.27^^ 
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2.3. Predictions for Counts and Redshift Distributions 

In Fig. 2, Fig.3 and Fig.4, predictions by these 3 models for the 15//m, 60/xm, 90//m, 
nOfim, SbO/im counts and for the CIB are compared with the data ^, respectively. The 
data points are taken from a large pool of measurements found in the literature, which is 
still evolving rapidly. This is particularly true for the measurements involving ISO data, 
because more and more new data reduction tools are becoming available for these rather 
complex data. 

Model SI, which is otherwise same as the 'peak model' in Xu et al. (2001) except for 
the modifications listed in Section 2.2, predicts that the contribution from the starburst 
population dominates almost everywhere in these plots. The model predictions are in 
reasonably good agreements with data in all plots, given the large dispersions in the data 
sets. Particularly, the new results of Gruppioni et al. (2002) on the 15/xm counts of ELAIS 
sources (filled 4-point stars in plot) are about a factor of 3 lower than the counts from 
the same survey reported by Serjeant et al. (2000) (open 4-point stars), and are also 
significantly lower than other ISOCAM measurements in the flux range of 1 mJy < fi^^m < 
10 mJy. Model SI predicts a less prominent sub-mJy peak in the 15fim counts than the 
'peak model' of Xu et al. (2001), because the sharp peak at z—1.5 in the evolution functions 
of the old model is replaced by a smooth peak in the new model (at z ~ 1.3, see Fig.l). It 
is worth to note that, in the bright end of the SCUBA counts {fsbOfim ~ 5m Jy), the model 
predicts that the contribution from galaxies with AGN exceeds that from starbursts. This 
result seems to be consistent with the hmited knowledge we have about the bright SCUBA 

^The units used in the plots throughout this paper are rather heterogeneous. This is 
becuase many data points collected from the literature are measured from figures that used 
different units. Often it is difficult to convert them to the same units, in particular for errors 
bars. 
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sources: out of the 7 SCUBA sources of fsbOum > SmJy in the sample of Smail et al. (2002), 
4 show signs of AGN. On the other hand, the model predictions on the SCUBA 850//m 
counts are close to the lower boundary of the data. 

In contrast to Model SI, Model S2 predicts that for 15/xm counts, 850/xm counts and 
for the CIB, the population of normal galaxies dominate (Fig.3). This is in agreement 
with the results of Rowan- Robinson (2001). The same model predicts that normal galaxies 
and starburst galaxies give about equal contributions to the 170//m counts, and starbursts 
dominate the 90/i.m and GO/^m counts. This model can also fit the data in all these plots 
well (Fig.3). 

By design, the contributions from normal galaxies and from starburst galaxies are 
indeed nearly equal in the predictions of Model S3 for 15/xm counts (Fig.4). The same 
model predicts that the dO/im, 90//m and 170//m counts are dominated by starbursts. For 
the SCUBA 850/xm band, the model predicts comparable contributions from the normal 
galaxies and from the starbursts to the counts in a wide flux range, and the contribution 
from AGNs dominant the counts at a few 10s of mJy level while being neghgible at 
sub-mJy level. The CIB predicted by this model is dominated by the contribution from the 
starbursts in the wavelength range 20/im ^ A ^ 300/i.m, and by that from normal galaxies 
at other wavelengths. 

In Fig. 5, Fig.6 and Fig. 7, predictions by the 3 models for the redshift distributions for 
the IRAS 60/xm, ISOCAM IS^im, ELAIS QO^im, FIRBACK 170//m, and SCUBA 850//m, 
and future SWIRE 24//m surveys are plotted. The model predictions for the IRAS 60//m 
and ISOCAM 15/xm surveys are compared to the data. All 3 models give more satisfactory 
fits to redshift distribution of the IRAS 60/xm sources than the 'peak model' (favorite 
model) of Xu et al. (2001) does, though Model SI still over-predicts the number of IRAS 
dOfjim sources with z > 0.2 by about a factor of 3. For the redshift distribution of ISOCAM 
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15/xm sources in the HDF-north, the predictions of the Model SI gives the best fit among 
the 3 models although, given the large error bars of the data, the difference among model 
predictions is subtle. 

All the 3 models predict bi-modal redshift distributions for ELAIS 90/xm and 
FIRBACK nOiim sources, with normal galaxies dominating the first peak and starburst 
galaxies dominating the second peak. Model SI predicts that most of ELAIS 90//m and 
FIRBACK 170//m sources have z > 0.5. In contrast. Model S2 predicts that most of these 
sources have z < 0.5, being in the first peak of the redshift distributions. This difference is 
due to the so-called ' temper ature-redshift' degeneracy for infrared galaxies (Blain 1999): 
In the long wavelength bands such as the 90//m and the 170//m, the sources in Model SI 
tend to have larger redshifts than those in Model S2 in order to compensate their warmer 
IR SEDs. Interestingly, this indicates that we can break the degeneracy over these models 
by comparing the predicted and observed redshift distributions of the ISOPHOT sources in 
these bands. Serjeant et al. (2001) obtained high confidence redshifts for 16 (out of 37) 
sources of /go^m > O.lJy in the ELAIS SI field (3.96 deg^), none of them has redshift larger 
than 0.5. This sets a lower limit of 43±11% for the percentage of the 90/xm sources in this 
field to have z < 0.5. Taken at the face value, only the prediction of S2 (53%) is comfortably 
above this lower limit, while the prediction of S3 (37%) is marginally consistent with it, and 
the prediction of SI (30%) is marginally below the limit. The above comparisons seem to 
favor Model S2. However, given the fact that this data set was taken from a small region 
and can be seriously affected by any clustering effect, more observations are needed for any 
definitive conclusions. 

Model SI predicts a broad redshift distribution, peaking between 2 < z < 3, for bright 

SCUBA sources (/sso^tm > 8 mJy), broadly consistent with observational constraints (Small 
et al. 2002; Ivison et al. 2002; Fox et al. 2002). For the same sources. Model S2 predicts a 
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prominent peak around z—1.5 and less than 25% of sources at redshifts > 2, not favored by 
observations which in general suggest a mean redshift larger than 2 (see, e.g. Ivison et al. 
2002). Model S3 also predicts a peak around z=1.5 for the /sso/im > 8 mJy sources, though 
with a much wider high- redshift wing (44% of sources have z >2). 

2.4. Predictions for ISO LFs and SFH 

The ISOCAM 15/xm sources in HDF-North have been thoroughly studied in the 
literature (Rowan- Robinson et al. 1997; Aussel et al. 1999; Elbaz et al. 1999; Cohen et al. 
2000; Franceschini et al. 2001; Elbaz et al. 2002). All but 1 of a total of 41 sources with 
fibum > 0.1 mJy have spectroscopic redshifts (Franceschini et al. 2001; Elbaz et al. 2002). 
Exploiting these redshifts, we derived the 15//m luminosity functions (LFs) in three redshift 
intervals: 0.4 < z < 0.7, 0.7 < z < 1, and 1 < z < 1.3 using the classical Vmax method. 
These results are compared to the predictions of the three models in Fig. 8. 

Given the small size of the data set (40 galaxies), it is understandable that the derived 
luminosity functions have substantial uncertainties. Particularly, significant fluctuations 
can be found in the redshift distributions (Fig.5-7), presumably due to clustering effects. 
Compared to these uncertainties, the difference in the predictions by the three different 
models is again subtle. They all lie near the lower ends of the error-bars (Poisson error) 
of the two data points in the panel of 0.4 < 2; < 0.7. This is likely to be due to a bias in 
the data caused by the over-density in the 0.4 < 2; < 0.6 bin (Fig.5-7). In the panels of 
0.7 < z < 1 and 1 < z < 1.3, the predictions of Model SI are slightly higher and fit the 
data slightly better than those of other two models. 

In Fig. 9a, the star formation history (SFH) predicted by the three models are compared 
to observations. The red lower limits are derived from the ISOCAM Ib/im LFs (Fig.8). 
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The lower limit of Ibjim luminosity density is derived by summing up the contributions in 
the bins where the LF is actually measured (i.e. no extrapolations are applied to measured 
LFs), then converted to total IR luminosity using the formula L^j. = 11.1 x I/is^^ (Elbaz et 
al. 2002). The star formation rate is then estimated from the IR luminosity density using 
the formula of Kennicutt (1998): SFR (M© yr-^)=Li^ x 1.7 10"^° L©. The predictions of 
the three models are very close to each other (within 50%). 

In Fig. 9b, Fig. 9c and Fig.9d, these predictions are broken into contributions by different 
populations. Indeed Model SI predicts that, except for in the low z end {z ^ 0.2), the 
starburst galaxies overwhelmingly dominate the star formation in the universe. Model S2 
predicts just the opposite: the star formation in the universe has been always dominated 
by normal late-type galaxies, while in Model S3 the contributions of these two populations 
are more comparable. All three models predict minor contributions from galaxies with 
AGNs. Note that since some of the IR emission from galaxies with AGNs is powered by 
the gravitational energy released in the AGN, not by star formation, the model predictions 
plotted here should be treated as upper limits. 

In the literature, the steep decline of the star formation rate since z ~ 1 has been well 
established since early works of Lilly et al. (1996) and Madau et al. (1996). However, 
it is still controversial on what happened to the star formation rate in earlier universe, 
particularly before z — 2. Most information on this issue is obtained from observations of 
Lyman-Break galaxies (LBGs, Steidel et al. 1999). However, the extinction correction for 
these rest-frame UV selected galaxies is very uncertain. In addition, many high redshift 
star-forming galaxies may be completely missed in the surveys of LBGs because of heavy 
extinction, causing systematic underestimation when using LBGs to determine the star 
formation rate in high z universe (see, e.g., Rowan-Robinson ct al. 1997). In this respect, 
the predictions of our models follow closely the data points determined from LBG surveys 
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after the extinction correction (Steidel et al. 1999). As discussed in Xu et al. (2001), the 
evolution of IR galaxies at redshifts > 2 is mainly constrained by the sub-mm data in 
the SCUBA surveys and in the CIB observations, due to the negative K-correction. Too 
much star formation in those large redshifts will result in too many sub-mm counts and 
too high sub-mm background radiation. In the literature, other models (e.g. Gispert et 
al. 2000; Rowan- Robinson 2000; Chary & Elbaz 2001) which also invoked sub-mm data 
to constrain the evolution of ^ > 2 IR galaxies found the similar trend (a peak at ^ ~ 1 
and a shallow/flat slope in z ^ 2) in the SFR evolution. The claim that the universal star 
formation rate increases monotonically with redshift to very high z ( ^ 7) by Lanzetta et 
al. (2002) is not supported by our and other studies on evolution of IR galaxies. 



3. Passive Evolution Models for E/SO Galaxies 

The models follow a simple, passive evolution approach (Pozzetti et al. 1996). The 
basic assumption is that there has been no star formation in an E/SO galaxy since its 
initial formation. Consequently, its radiation in different bands (i.e. the SED and the L/M 
ratio) evolves passively with the ever-aging stellar population. Instead of assuming that all 
E/SO's formed at once together (as in the classical monolithic galaxy formation scenario), 
the E/SO galaxies are assumed to form in a broad redshift range (Pranceschini et al. 1998), 
specified by a truncated Gaussian function. The SEDs of different ages arc calculated using 
the GRASIL code of Silva et al. (1998). No dust emission is considered in these models. 



3.1. Evolution Function 



The evolution function ^(M, i), specifying how many E/SOs are formed in a unit 
volume, in a unit time interval, and in a unit mass interval, is assumed to have the following 
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form: 

^(M,t) = ^(M) T(t) (4) 
where the time dependence function T{t) is a truncated Gaussian: 

T{t) = exp{-{t-t{Zpeak))V0U^) {t > t{zo)) 

= (t<t{zo)). (5) 

In this prescription, for the sake of simphcity, the time dependence is assumed to be 
independent of the mass. This imphes that E/SOs of different mass have the same formation 
history (i.e. no differential evolution). The time dependence function T{t) is fully defined 
by 3 free parameters: the peak formation time t{zpeak), the time scale lu (in Gyr), and the 
time when the formation of E/SOs started t{zo). 

Accordingly, for a given redshift z, the number of E/SO's in a unit mass interval and a 
unit volume can be found from the following integration: 

$(M, z) = ip{M) / T{t') dt'. (6) 

Jt{zo) 

If the local mass function $(M, 2; = 0) is known, then the mass dependence function 
ip{M) can be derived: 

^(M) = . (7) 

where t{z = 0) is the current age of the universe. 

In this work, the local mass function is constrained using the local K-band luminosity 
function of Kochanek et al. (2001) and the Lk/M ratio predicted by GRASIL code (Silva 
et al. 1999). 
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3.2. Galctxy Age Distribution and SED Assignment 

Eq(6) predicts galaxy number for a given redshift in a unit mass interval and a unit 
volume. These galaxies have different ages spanning between t — and r = t{z) — t{zo). 
The age distribution function, G{z,t), is also a truncated Gaussian: 

G{z,t) = Goexp{-{t{z)-T-t{zj>,ak)f/^'') 0<T<{t{z)-t{zo)) 

— (otherwise). (8) 

For every galaxy simulated, an age is assigned to it according to the above distribution 
(i.e. probability) function. According to the passive evolution model, galaxies of different 
ages have different SEDs and different L/M ratios. The SEDs of different ages are calculated 
using GRASIL (Silva et al. 1998). Again, for the sake of simplicity, we assume a simple 
'single-burst' scenario, in analogy to the merger scenario of E/SO formation as hinted 
at in the studies of ULIRGs (Kormendy & Sanders 1992), to model the formation of 
E/SOs. Namely, we assume that all the stars in an E/SO galaxy were formed in a short 
burst (lasting 10® yrs), after which all the ISM was blown out. Accordingly, the following 
parameters are adopted for the input of GRASIL: t^in = 0.1 (Gyr), ksch = 1-0, I'sch = 20.0, 
and Tinf — 0.01. Then, in the simulation of E/SOs, we include only sources older than 1 
Gyrs. Here we implicitly assume that merger remnants younger than 1 Gyr do not look like 
E/SOs because these systems may have not fully relaxed. Namely they should be classified 
as post-starbursts or E+As, not E/SOs. Optical follow-up observations of ISOCAM sources 
(Flores et al. 1999b; Aussel et al. 1999; Cohen et al. 2000; Elbaz et al. 2002) have shown 
that E+A galaxies are IR bright and have similar IR properties as active starburst galaxies. 
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3.3. Three Models for E/SO Evolution 

We consider 3 different E/SO models in this paper, specified by different Zpeak and 
uj (Table 2, Fig. 10). Model El {zpeak = 5, a; = 0.5 Gyr) mimics the classical monolithic 
scenario (Eggen et al. 1962). Model E2 assumes a later and broader E/SO formation epoch 
{zpeak — 2, Lo — 2 Gyr). Model E3 {zpeak — 1, CO — 3 Gyr) assumes that most of E/SOs are 
formed at z > 1, as hierarchical galaxy formation workers have advocated (Kauffmann et 
al. 1993; Kauffmann & Chariot 1998). 

Table 2. Evolution models for E/SOs 

Models Zpeak (Gyr) zq 

El 5 0.5 7 

E2 2 2 7 

E3 13 7 

In Fig. 11 we compare the model predictions on the redshift dependence of optical/NIR 
colors with the data from Franceschini et al. (1998). Predictions of all three models can 
fit the overall trend of the data well. On the other hand, Model El predicts very little 
scattering, which is not confirmed by the data. The other two models predict significantly 
larger dispersion, in better agreement with the data. 

In Fig. 12 the redshift distribution of early type galaxies of Kab < 20.15, taken from 
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Pranceschini et al. (1998), is compared to the model predictions. Note that only 15 out 
of 35 redshifts in Pranceschini et al. (1998) are spectroscopically measured, the rest are 
photometric redshifts. Among the three models. Model E2 fits the data best. It should be 
noted that, as discussed in Benitez et al. (1999) and Rodighiero et al. (2001), some high z 
E/SOs may have been missed in the Pranceschini et al. (1998) investigation. 

4. Counts in UV, Optical, NIR and MIR Bands 

Both the E/SO galaxies and the dusty galaxies contribute significantly in these bands. 
Since the three E/SO models predict almost identical counts (difference ^ O.l dex), in the 
rest of the paper we shall present only the results from Model E2 {zpeak — 2, u; — 2 Gyr). 

In Pig. 13, Pig. 14 and Pig. 15 the contributions of E/SO galaxies to the counts in these 
bands (as predicted by Model E2) are added to the contributions of dusty galaxies as 
predicted by Model SI, Model S2 and Model S3, respectively. These model predictions 
are compared to observations in the vacuum UV (2000A), the B and R, the NIR K, and 
the MIR 6.7/xm band (ISOCAM) and 12^m (IRAS and ISO CAM) bands. Between the 
three figures, the differences in the total counts in the bands plotted are generally within 
the uncertainties of data. The model predictions can account for all counts in the R, K 
and MIR bands. On the other hand, the model predictions are significantly lower than 
the counts observed in the 2000A band (Milliard et al. 1994). Por the B band, the model 
predictions are shghtly lower (~ 0.3 dex) than the observed counts for S > 20 mags. These 
results strongly hint at a population of star forming galaxies which are infrared quiet (very 
low dust attenuation/emission), therefore seen only in the UV and the blue bands (see 
Section 7.3 for more discussion). 

In the bright end of the K-band counts, the prediction by S1+E2 (Pig. 13) is in good 
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agreement with the morphologically segregated counts reported by Huang et al. (1998), 
namely at K ^ 17 the E/SO galaxies contribute about 50% of the counts, and at fainter 
magnitudes the contribution from late- types becomes more and more dominant. The 
predictions by S2-I-E2 and S3-I-E2 on the E/SO contribution to K-band counts at K= 16 are 
~ 30% and ~ 40%, shghtly less than the observational result of Huang et al. (1998). 

The three models fit the faint, sub-mJy ISOCAM 6.7/xm counts very well. However, 
predictions of S1-I-E2 are about a factor of 3 lower than the ELAIS counts (Serjeant et al. 
2000) at few mJy level. For these counts, predictions by S3+E3 have the best agreement 
with, though being still about 50% less than, the data. It should be noted that, as indicated 
by the large discrepancy between the ELAIS 15//m counts of Serjeant at al. (2000) and of 
Gruppioni et al. (2002), the uncertainties of the ELAIS 6.7 iim counts of Serjeant et al. 
(2000) may be signicantly larger than reported. 

The scatters of the ISOCAM 12//m data (Clements et al. 1999) are large, indicating 
significantly larger uncertainties than the Poisson noise (error bars in the plot). The three 
brightest ISO points may suffer serious biases caused by errors in galaxy/star separation 
(Clements et al. 1999), so are hkely to be less rehable than other data points. If this is 
indeed the case then, among the three models, the predictions of S1+E2 have the best 
agreement with the data. The effect of the Local Supercluster (Lonsdale et al. 1990) can 
be seen in the bright IRAS 12//m counts. 

5. Counts and Confusion Limits in SIRTF Bands 

The E/SO galaxies contribute significantly only to the four IRAC bands (3.6//m, 4.5//m, 
5.8/im, and 8/i.m), while the 3 MIPS bands (24/i,m, 70/im and 160/im) will see only dusty 
galaxies. 
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Model predictions of contributions by E/SO galaxies and by different populations of 
dusty galaxies to the counts in three SIRTF bands (3.6/im, 24/xm and bands) are 

plotted in Fig. 16. The total counts by the three different models are very close to each 
other, though the relative contributions from normal galaxies and from starbursts are very 
different. This can be understood by the fact that, tuned to fit the same counts in the ISO 
bands which covers a wide wavelength range from to 170//.m, these models are forced 

to be similar with each other. 

The 3-0" confusion limits are given in Table 3 for all 7 SIRTF bands. These are 
calculated by the method described in Xu et al. (2001), and assuming that SIRTF (85 
cm dish) is diffraction limited in all bands (so the beams can be approximated by the 
Airy function). Since this idealized assumption may not be true, particularly for the short 
wavelength IRAC bands (e.g. the 3.6/xm and 4.5/xm bands), the results for those bands 
should be treated as lower hmits. The confusion hmits predicted by the 3 different models 
differ by up to 80%. The largest difference occurs for the 70//m band. These differences 
reflect the real uncertainties, mostly due to the uncertainties in the ISO data which are the 
major constraints to the models. 



Table 3. Confusion limits {3a) of SIRTF bands 

Models 3.6/xm i.biim 5.8/im 8/im 2Aiim IQiim 160/xm 

S1+E2 0.12 //Jy 0.21 //Jy 0.43 //Jy 1.02 //Jy 65 /xJy 6.15 mJy 81.5 mJy 

S2+E2 0.18 //Jy 0.35 //Jy 0.65 //Jy 1.32 //Jy 57 //Jy 3.35 mJy 57.5 mJy 

S3+E2 0.14 /xJy 0.27 /xJy 0.53 /xJy 1.17 /xJy 57 /xJy 4.45 mJy 68.5 mJy 
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6. Colors as Model Discriminators for Dusty Galetxies 

In empirical evolution models, different populations of dusty galaxies are distinguished 
by different SEDs. Therefore, the most direct method to distinguish models predicting 
different dominant populations for IR sources is to compare the predicted color distributions 
with the observations. 

In Fig. 17, we compare predictions of the three models for evolution of dusty galaxies 
(Table 1) for colors of ISO galaxies. The detection limits are set to be /e.r^m > 30 // Jy, 
/i5/im > 100 nJy, fgofxm > 100 mJy, firoixm > 180 mJy, i? < 24 mag, and K < 20 mag. 
The predictions for /i5^m//i? and fi^^m/ Ik are compared to the data of ISOCAM 15/i.m 
sources in HDF-N (Aussel et al. 1999; Cohen et al. 2000; Hogg et al. 2000). It appears 
that, for ISO galaxies, only the fisi^m/ fx color is a good model discriminator in the sense 
that the peaks of the color distributions predicted by the three models are separated from 
each other. The data seem to favor Model SI which predicts a peak in the fi^^m/ Ik color 
distribution right at the place where the data peaks. However, this data set is again too 
small (30 galaxies) and error bars too large to distinguish the three models. 

SWIRE^ a SIRTF Legacy Science program, will survey 65 deg^ of sky in all 7 SIRTF 
bands (3.6, 4.5, 5.8, 8, 24, 70, and 160//m) down to a depth comparable to, or even deeper 
than, that of ISOCAM 15/im surveys. The documented 5a sensitivity limits of SWIRE are: 
h.&nm = 7.3 nJy, U^^rn = 8.7 iiJy, /s.g^,^ = 27.5 /iJy, /g^^ = 32.5 /iJy, /24^„ = 0.45 mJy, 
and Iroixm — 2.75 mJy, /leo/im = 17.5 mJy. Compared to predicted confusion limits in 
Table 3, it appears that SWIRE surveys will be confusion limited in the the 70/xm and 
IGO/im bands. Extensive ground based follow-up observations in the optical, NIR and radio 
bands will be carried out. These survey areas will also be observed in the far-UV (1500A 



^http://www.ipac.caltech.edu/SWIRE. 
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and 2300A) by GALEX^ in the deep survey mode. In Fig. 18, model predictions for the 
distributions of 6 colors of SWIRE galaxies are plotted. These colors are selected among 
many possible combinations to illustrate how the colors of SWIRE galaxies can discriminate 
the models. Model simulations of sky coverage of 5 deg^ are carried out. For each plot, the 
samples of simulated sources are selected according to SWIRE's sensitivity hmits or, for 
the 70/um and 160/im bands, the confusion limits. Also it is required that R < 24 mag and 
K <20 mag when the R and K band data are invovled. In four {froixm/ fs/xm, f24^im/ fssum, 
/24/im//8/um, ^ud f24nm/ fx) of the six colors plotted, the peaks of Model SI (starburst 
dominant) and Model S2 (normals dominant) are clearly separate. Among these, the peak 
of the distribution predicted by Model 3 (intermediate model) is close to that of SI in the 
fiOfiml fsfim color plot, closc to that of S2 in the fiAixm/ fsfim color plot, and more ambiguous 
in the other two plots. Given the large sky coverage of SWIRE, which is 13 times of what 
is simulated here, it is very hopeful that these color distributions will indeed provide clues 
to the question of which population is dominant among IR sources. 

7. Discussion 

7.1. Evolution of SEDs of Steirburst Galcixies 

As pointed out in Xu et al. (2001), the most important assumption in our models 
for the evolution of dusty galaxies is that high redshift (i.e. z ^ 1) star-forming galaxies 
have the same SEDs as their local counterparts when the luminosity is the same. Since 
our SEDs cover from the UV all the way through to the radio waveband, the validity 
of this assumption demands similarity between high redshift and local galaxies in many 
physical conditions, a requirement that seems too strong to be fulfilled in the strict sense. 



^http: / / www.srl.caltech.edu/galex. 
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Particularly, the long term star formation history plays an important role in the optical and 
NIR emission of galaxies, and it is obvious that high redshift galaxies have very different 
long term star formation history (i.e. much younger) than that of local galaxies. Does this 
mean that our predictions for the optical/NIR flux densities of high redshift dusty galaxies 
are flawed and therefore unreliable? Our argument to dispute this suspicion is based on 
the fact that high redshift galaxies (particularly IR selected galaxies), already detected 
or to be detected in the future surveys, are almost exclusively high luminousity galaxies 
(less luminous galaxies with high redshift are too faint to be detected). As their local 
counterparts, these high redshift ULIRGs must host very powerful starbursts or AGNs as 
energy sources which contribute much of the emission even in the NIR bands (~ 40% for the 
local ULIRGs, Surace et al. 2000; Scoville et al. 2000). Therefore, as far as the detectable 
IR sources are concerned, the difference in the underlying old population between the 
high redshift dusty starforming galaxies and their local counterparts does not affect very 
seriously our model predictions. 

Recent multi-wavebands observations of SCUBA galaxies (Small et al. 2002), 
Lyman-Break Galaxies (Adelberger & Steidel 2000), and ISOCAM galaxies (Flores et 
al. 1999a; Aussel et al. 1999; Cohen et al. 2000; Elbaz et al. 2002) are consistent 
with our assumption that high redshift star-forming galaxies have similar SEDs as their 
local counterparts, particularly when binned according to the luminosity (as stressed by 
Adelberger & Steidel 2000). On the other hand, some new observations hint at possible 
systematic differences between high redshift and local ULIRGs. The HST image of SCUBA 
galaxy SMM J 140 11 +0252 (Ivison ct al. 2001) shows that the star formation activity in 
that source is widely spread (up to a few kpc) , a situation remarkably different from typical 
ULIRGs found in the local universe (Sanders & Mirabel 1998). Compared to centrally 
concentrated starbursts found in most of the local ULIRGs, galaxies with such widely 
distributed starbursts are expected to have less steep MIR slopes (i.e. smaller f25nm/ fi2nm 



-26- 



ratio) and cooler FIR color (i.e. larger fwoi^m/ feoum ratios). This is due to the less intense 
radiation field (less warm dust emission in the 25/xm and 60/^m bands) and less dust 
opacity (less extinction for MIR fluxes). Indeed, Chapman et al. (2002) found that two 
sources detected both by the FIRBACK 170//m band survey (Puget et al. 1999; Dole et 
al. 2001) and by SCUBA (Scott et al. 2000), one at z=0.91 and the other at z=0.46, have 
significantly cooler dust temperatures {T^ust ~ 30K) compared to T^ust ~ 50K found for 
typical ULIRGs such as Arp220. They argue that this may indicate the starbursts in these 
systems are also extended. Given the small amount of information and the possible bias for 
cooler galaxies due to the sub- mm selection (Chapman et al. 2002), it is still too early to 
tell whether high-redshift ULIRGa are systematically more extended, and therefore have 
SEDs closer to less luminous local interacting galaxies (such as the Antennae galaxies) 
which are in earher stages along the merger sequence. Future surveys hke SWIRE will 
address these questions. 



7.2. E/SO evolution and starburst evolution 

Comparisons between predictions E/SO evolution models and observational data 
(Fig. 11-12) indicate that the intermediate model (Model E2), which assumes a peak 
formation redshift of Zp^ak = 2, is the favorite among the three models. This is in agreement 
with some previous works (Franceschini et al. 1998; Rodighiero et al. 2001; Im et al. 
2002). This model can reproduce well the trend in the color-z plots and also fit the 
redshift distribution very well. However, particular in the B-K v.s. z plot, the data show 
significantly larger dispersion than the model predictions. There are two possible causes for 
this: 

(1) In our simple models, we have assumed that all of the E/SOs are formed through 
the same starburst procedure, which means their stellar populations have the same 
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metallicity. It is known (e.g. Worthy 1994) that local E/SOs have different metaUicity 
which is a major cause of the different M/L ratios and colors among these galaxies. 
Therefore, by neglecting these effects, our models leave some of the scatters in the 
color distributions unaccounted for. 

(2) Some of the blue E/SOs in the data (Franceschini et al. 1998) may not be true E/SOs. 
According to Im et al. (2002), many of these 'blue interlopers' have strong, narrow 
emission lines, suggesting that they are low-mass starbursts rather than massive 
star- forming E/SOs. 

In fact, there is still a lack of consensus in the definition of E/SOs in deep surveys. 
Using a strict algorithm selecting the most symmetric and smooth galaxies, Im et al. 
(2002) found far less blue sources among their E/SO sample than, e.g., Schade et al. (1999) 
and Menanteau (1999) whose samples were selected with less strict algorithms. A related 
uncertainty in our models is the choice on the exclusion of galaxies younger than 1 Gyr 
(Section 3.2). This choice is not entirely arbitrary: pushing the cut-off toward younger 
ages means more galaxies with high L/M ratios, which in turn results in too high model 
predictions for optical and NIR counts. This shows that the question of where to put the 
boundary between E/SOs and post-starbursts deserves more investigation. 

If indeed E/SOs are formed through mergers (Toomre 1978; Kormendy & Sanders 
1992), then their evolution is linked to the evolution of star burst galaxies which are 
closely related to mergers (Sanders & Mirabel 1998). Our results indeed indicate some 
synchronization between the two populations, in the sense that the peak of the formation 
function of E/SOs in the best fit model {zpeak — 2) is close to the peak of the evolution 
functions of starburst galaxies {zpeak — 1-4). Future works exploring this possible link will 
provide important constraints to the evolution of both populations. 



-28- 



7.3. IR-quiet stair-forming galatxies 

Our results (Fig. 13 - 15) show that only a small fraction (10 - 20 %) of UV selected 
galaxies (as in the sample of FOCA survey, Milliard et al. 1992) are IR bright, as indicated 
by the small contribution of simulated dusty galaxies to the UV counts (the E/SO galaxies 
contribute even less). This suggests a separate population of IR quiet, UV bright galaxies 
which dominate the UV selected samples. Such a population is also needed for the B-band 
counts, because our simulations under-predict 30 - 50% of observed counts in that band, 
while fully account for the R and K band counts (Fig. 13 - 15). The best candidates 
for such galaxies are the low-metallicity (therefore low dust content) blue dwarf galaxies 
such as I Zw 18 (Searle & Sargent 1972). There is an apparent link between this IR quiet 
star-forming galaxy population and the 'faint blue galaxies' found in deep optical surveys 
(see Koo & Kron 1992 and EUis 1997 for reviews). What is the relation between this 
population and the dusty (IR bright) star-forming galaxies? How does this population 
evolve (backwardly) with the redshift, and how does this evolution correlate with the 
evolution of IR bright galaxies? Are Lyman Break Galaxies, being selected by the UV flux 
in the rest frame, more closely related to the IR quiet star forming galaxies, or to the dusty 
star forming galaxies (as argued by Adelberger & Steidel 2000)? Answers to these questions 
will help to unify the pictures of galaxy formation/evolution seen in different wavebands. 
We plan to address these questions in our future work, particularly in connections with the 
forthcoming GALEX and SWIRE missions. 

8. Summciry 

New models for the evolution of extragalactic IR sources are presented in this paper. 
The models for dusty galaxies and for E/SO galaxies, the latter contributing significantly to 
counts at wavelengths A < 10/xm, have been developed separately. 
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Compared to previously published models in Xu et al. (2001), the new models for 
evolution of dusty galaxies in this work have the following improvements: 

(1) The evolution functions have the form of smoothly-joined 3-piece power-law (Fig.l), 

instead of the sharply- joined 2-piece power-law. 

(2) New local luminosity functions at 25/xm, which take into account the evolution 
effects, are used for the three dust galaxy populations (i.e. normal late-type galaxies, 
starburst galaxies, and galaxies with AGNs). 

(3) The UV portion (A < SOOOA) of the SEDs in the SED library is constrained by an 

empirical correlation between fuv — Jb and B — K, instead of mere an extrapolation 
of the optical SED (A > 4000A). 

In order to address the question whether normal late-type galaxies or starburst galaxies 
dominate among IR sources oi z > 0.5, three new models are developed: 

• Model SI — starburst galaxies dominant ; 

• Model S2 — normal galaxies dominant; 

• Model S3 — intermediate between SI and S2. 



Predictions of these three models for counts in various bands are fairly close to each other, 
therefore they can hardly be distinguished using the counts. They can also fit very well 
the luminosity functions of ISO IB/im sources in three redshift intervals (0.4 < z < 0.7, 
0.7<2;<1, 1<2;<1.3). In principle, they can be distinguished by redshift distributions 
of FIR sources, but the large beams of FIR detectors (such as the ISOPHOT cameras) 
make it rather difficult to pin-down the optical counterparts (usually very faint) of faint 
FIR sources. At the same time, the peaks in the distributions of several IR and optical 
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colors predicted by these models have separated locations. We argue that these color 
distributions are the best tools to distinguish these models. There are only very limited 
amount of multi-waveband data available for high redshift dusty galaxies in the literature 
(mostly for ISOCAM 15/im sources), a situation that will be drastically improved when 
SIRTF is launched and SWIRE surveys are available. 

The models for E/SOs follow a simple, passive evolution approach. The basic 
assumption is that there has been no star formation in an E/SO galaxy since its initial 
formation. Instead of assuming that all E/SO's formed at once together (as in the classical 
monolithic galaxy formation scenario), the E/SO galaxies are assumed to form in a broad 
redshift range (Pranceschini et al. 1998), specified by a truncated Gaussian function. 
The SEDs of different ages are calculated using the GRASIL code of Silva et al. (1998). 
No dust emission is considered in these models. Three such models with different E/SO 
formation histories are calculated: Model El {zpeak — 5, u — 0.5 Gyr) is close to the 
classical monohthic scenario. Model E2 assumes a later and broader E/SO formation epoch 
{zpeak = 2, uj = 2 Gyr). Model E3 {zpeak = 1, = 3 Gyr) assumes that most of E/SOs are 
formed at z > 1, as the hierarchical galaxy formation works have advocated (Kauffmann 
et al. 1993; Kauffmann & Chariot 1998). Comparisons with limited data (e.g. colors 
and redshift distribution) available for morphologically identified E/SO galaxies at z~ 1 
indicate that model E2 can fit the data best among the three models. This suggests a 
synchronization between the evolution of E/SO galaxies and of starburst galaxies, in the 
sense that the peak of the formation function of E/SOs {zpeak — 2) is close to the peak of 
the evolution functions of starburst galaxies {zpeak = 1-4). Combining model predictions by 
E2 with those by SI, S2 and S3 (dusty galaxy evolution models), comparisons with number 
counts in different wavebands (Fig. 13-15) indicate that E/SOs contribute upto 30 - 50% of 
the optical/NIR counts in the bright end, and about 20 - 30% of the ISOCAM 6.7/xm band 
counts. Their contributions to counts in the UV (2000A) and in the longer wavelength IR 
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(> 12 jim) bands are negligible. 

Using these new models for extragalactic IR sources, particularly including E/SO 
galaxies, we made new predictions for the counts and confusion limits in the SIRTF bands. 
The results indicate that SWIRE surveys will be confusion limited in the 70/xm and 160/xm 
bands. The confusion limits predicted by different models differ by up to 80%. These 
differences reflect the uncertainties of these predictions. 
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Appendix 

A. New 25/xm Local Luminosity Functions for Three Populations 

In Xu et al. (2001), luminosity functions were presented for three population 
subsamplcs-AGNs, normal late-type galaxies and starbursts. The subsamplcs were 
determined by IRAS colors. We have extended this work to take into account the effects 
of luminosity evolution on the LLF estimates, and to include a more comprehensive 
uncertainty analysis. 

We recomputed the population LLPs including a luminosity evolution term of the 
form (1 + zY, for values of q =0, 3.0, and 4.5. This factor is applied to both the source 
luminosity, and the minimum luminosity detectable at the source's redshift. The shape 
parameters are tabulated in Table A.l. 
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Table A.l. Pciremeters of 2hij,m Luminosity Functions 



Population 


q 


a 


/5 




a + j3 


AGNs 


0.0 


0.336 


1.691 


6.9x10*^ 


2.027 


AGNs 


3.0 


0.329 


1.713 


6.6x10^^ 


2.042 


AGNs 


4.5 


0.326 


1.724 


6.5x10^ 


2.040 


Starbursts 


0.0 


0.265 


2.283 


7.9x10*^ 


2.548 


Starbursts 


3.0 


0.265 


2.275 


7.7x10^^ 


2.540 


Starbursts 


4.5 


0.264 


2.300 


7.9x10^ 


2.564 


Normals 


0.0 


0.482 


3.875 


5.7x10^^ 


4.357 


Normals 


3.0 


0.480 


3.992 


5.8x10^ 


4.472 


Normals 


4.5 


0.479 


4.055 


5.8x10*^ 


4.534 



The resulting LLPs (with same normalization constant applied) are plotted as 
visibilities in panel a of Pig. 19 (sohd hne is no evolution; dashed is q=3.0; dashed-dot is 
q=4.5). There is a small difference between the LLPs at large luminosities. The model 
calculations in this paper use the LLPs with q=3.0. 

We performed an analysis of the covariance of the fitted parameters, using the 
information matrix (Efstathiou et al. 1988). In panels 6, c, and dof Pig. 19, we have plotted 
the 68% confidence intervals in pairs of the parameters, as ellipses under the assumption of 
normally distributed uncertainties (Avni 1976; Press et al. 1992). In general, the parameters 
are not strongly correlated, except for a and j3 for AGNs. This correlation is explained by 
the distribution of this population at higher relative redshifts; the high-luminosity slope of 
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the luminosity function depends on the sum of a and (3, so the parameters are degenerate 
when a is not well determined at low luminosities. For the same reason, the confidence 
intervals for AGNs for a and are larger for the other populations. The variations in 
the parameters with evolution exponent are for the most part small compared to the 
uncertainties in the parameters. 

B. Extrapolation of SEDs to the UV Bands 

In Xu et al. (2001), the UV (1000 — 4000A) SEDs are extrapolations from data points 
in the B, J, H and K bands, and therefore are not well constrained. In this work, this is 
improved by introducing the following constraints: 

(1) The UV-B v.s. B-K correlation. The correlation is estabhshed using a sample of 
galaxies detected both in the vacuum UV bands (1500 A - 2500 A) and in the FIR 
bands (IRAS). The sample is taken from Xu & Buat (1995). The K-band magnitudes 
are found in the 2MASS database. In Fig.20, the UV-B v.s. B-K color-color 
plot for this sample is presented. The sample is divided into galaxies with AGNs 
(/60/xm//25/xm < 0.4), galaxies with IR excess (/60/xm//25/xm > 0.4 and L/j^/Ls > 0.5), 
and normal late-types {feoixm/ f25fj,m > 0.4 and Lfir/Ls < 0.5). Data of the three 
ULIRGs observed by Trentham et al. (1999) in UV using HST are also plotted. They 
follow the same trend of IR excess galaxies, though with larger scatters. We found 
that the trend for normal late-types can be well fitted by the following function 

log(i//^(2000l)/i//^ (44001) = -0.2 - 0.3 x {{B - K) - 2f {{B - K) > 2) 

= -0.2 - 0.1 X {{B -K)-2) {{B -K)< 2),(B1) 
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and for galaxies with AGNs and galaxies with IR excess, it can be well fitted by a 
two-step linear function 

log(i//^(2000l)/i//^(4400l) = -0.4 x ((B - K) - 2.) {{B - K) > 2) 

= {{B-K)<2). (B2) 

(2) The UV slope. We use the following relation between the UV slope (3 {Fx cx A'', 

1200A < A < 2600 A) and the L fir /Lb ratio, found by Calzetti et al. (1995), to 
constrain the slope of the UV SEDs between 1200A and 2600A: 

P = 1.12 X log(L/,,/LB) - 0.94. (B3) 

(3) Lyman-break. A sharp drop-off is imposed to the SEDs shortward of 912A: The flux 

density decreases a factor of 30 from 9 12 A to 700 A. There is a less steep drop-off, 
about a factor of 5, from 1200A to 912A, which is to mimic the effect of the Lya 
absorption. 

As an example, the model SED of M82 is compared to data in Fig.21. 

This research has made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. This work has made use 
of data products from the Two Micron All Sky Survey (2MASS), which is a joint project of 
the University of Massachusetts and the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronautics and Space Administration and 
the National Science Foundation. C. K. Xu, C. J. Lonsdale and D. L. Shupe were supported 
by the SIRTF Legacy Science Program provided by NASA through a contract with the Jet 
Propulsion Laboratory, California Institute of Technology. 
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Figure Captions: 



Fig. 1. — Schematic plots of evolution functions of the three new models. 
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Fig. 2. Comparisons of predictions of Model SI with observed IR counts and CIB. Data points in the plot of 15/im counts: 

the ELAIS counts reported by Serjeant et al. (2000) are plotted with open four-point stars, and the ELAIS counts reported 
by Gruppioni et al. (2002) arc plotted with filled four-point stars. Other data points have the same symbols as in Elbaz et al. 

(1999) ; A2390 (six-point stars); ISOHDF-North (open circles), ISOHDF-South (filled circles), Marano FIRBACK Ultra-Deep 
(open squares), Marano Ultra-Deep (exes), Marano FIRBACK Deep (asterisks), Lockman Deep (open triangles), Lockman 
Shallow (filled triangles). Filled squares with error bars are counts taken from Xu (2000). Note that the high points in the 
bright end (/is^m > 0.5 Jy) are due to Local Supercluster (Lonsdale et al. 1990). The shaded area marks the range of counts 
estimated by Maazei et al (2001). Data points in the IRAS 60/im plot: Large filled circles: Mazzei et al. (2001); Xs: Hacking & 
Houck (1987); open stars: Gregorich et al. (1995); open circles: Bertin et al. (1997); small filled squares: counts in the South 
Galactic cap {b" < -50°) by Lonsdale et al. (1990); open triangles: Saunders et al. (1991); open squares: Rowan-Robinson 
et al. (1990). Data points in the 90//m plot: filled circles: Efstathiou et al. (2000), crosses: Linden- Voernle et al. (2000), 
open square: Matsuhara et al. (2000), open diamonds: total counts of Juvela et al. (2000), open triangles: counts of multiple 
detections of Juvela et al. (2000). Data points in the 170/Ltm plot: filled circles: Dole et al. (2001), open square: Matsuhara et 
al. (2000), open diamonds: total counts of Juvela et al. (2000), open triangles: counts of multiple detections of Juvela et al. 

(2000) . The SCUBA 850;um counts: crosses: Blain et al. (1999); open circles: Hughes et al. (1998); open diamonds: Eales et 
al. 2000; open squares: Barger et al. (1999); filled diamonds: Scott et al. (2002). The cosmic IR background: Filled circles: 
Lagache et al. (1998); open squares: Finkbeiner et al. (2000); open stars: Gorjian et al. (2000); filled star: Dwek and Arendt 
(1998); large crosses; SCUBA source count results (Blain et al. 1999); shadowed area; the range of COBE/FIRAS results 
(Fixsen et al. 1998); diamonds and Xs with upper-limits: upper-limits from TeV gamma-ray radiation of Mrk403 and Mrk501 
(Dwek & Slavin 1994; Stanev & Franceschini 1998). 

Fig. 3. — Comparisons of predictions of Model S2 with observed IR counts and CIB. The 
sources of data points are the same as in Fig. 2. 

Fig. 4. — Comparisons of predictions of Model S3 with observed IR counts and CIB. The 
sources of data points are the same as in Fig. 2. 
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Fig. 5. — Predictions of Model SI for redshift distributions of various IR surveys. The 
observed redshift distributions (histograms in corresponding plots) of IRAS GO/im sources 
and ISOCAM 15//m sources are taken from Rowan- Robinson (2001) and Pranceschini et al. 
(2001), respectively. 

Fig. 6. — Predictions of Model S2 for redshift distributions of various IR surveys. Otherwise 
same as in Fig. 5. 

Fig. 7. — Predictions of Model S3 for redshift distributions of various IR surveys. Otherwise 
same as in Fig. 5. 

Fig. 8. — Comparisons of model predictions for the 15/xm luminosity functions (LFs) with 
observations. The data points for z—0 LF are taken from Xu (2000). The sohd hue is the 
Schechter function fit of the data points. Points in other panels are derived using data taken 
from Elbaz et al. (2002). 

Fig. 9. — Star formation history (SFH) predicted by new models. Data points: Red arrows 
with error bars: lower limits derived from the Ib/im LFs in Fig. 8 (see text). Filled squares: 
SFR from ISOCAM surveys (Flores et al. 1999b). Xs: SFR from Balmer line surveys 
(Gallego et al. 1995; Tresse & Maddox 1998; Yan et al. 1999; Glazebrook et al. 1999). Open 
squares: SFR from UV surveys ('extinction corrected' data taken from Fig. 9 of Steidel et al. 
1999). Open triangle: SFR from SCUBA (Barger et al. 1999). The shaded area: predictions 
by Pei & Fall (1995). All data have been converted to the cosmology model specified by 
Ho=75 km/sec/Mpc, flm = 0.3, and fl^ = 0.7. 

Fig. 10. — Schematic plots of formation functions of three E/SO models. 

Fig. 11. — Color v.s. redshift diagrams of E/SO galaxies. Model predictions compared with 
observations (open diamonds, Franceschini et al. 1998). 
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Fig. 12. — Model predictions for the redshift distribution of E/SOs in HDF-N compared with 
the data (histogram, Franceschini et al. 1998). 

Fig. 13.— Counts in the UV (2000A), B, R and K bands. Predictions of Model SI (for dusty 
galaxies) plus predictions of Model E2 (for E/SOs) are compared to data. Data points in the 
UV (2000A) plot: MiUiard et al. (1994). Data points in the B-band plot: Open squares: 
WiUiams et al. (1996); filled squares: Metcalfe et al. (1995); Xs: Metcalfe et al. (1991); 
open diamonds: Gardner ct al. (1996). Data points in the R-band plot: open squares: Lin et 
al. (1999); open diamonds: Cohen (2002). Data points in the K-band plot: open triangles: 
Bershady et al. (1998); open squares: Soifer et al. (1994); Xs: Minezaki et al. (1998); 
open diamonds: Gardner et al. (1996). Data points in the 6.7/xm band plot: open squares: 
Taniguchi et al. 1997; Xs: Altieri et al. 1999; diamond: Flores et al. 1999a; shaded area: 
Oliver et al. (2002) and Serjeant et al. (2000). Data points in the 12//m band plot: open 
squares: Fan et al. 1998; Xs: Clements et al. 1999. 

Fig. 14.— Counts in the UV (2000A), B, R and K bands. Predictions of Model S2 (for dusty 
galaxies) plus predictions of Model E2 (for E/SOs) are compared to data. Data points are 
the same as in Fig. 13. 

Fig. 15.— Counts in the UV (2000A), B, R and K bands. Predictions of Model S3 (for dusty 
galaxies) plus predictions of Model E2 (for E/SOs) are compared to data. Data points are 
the same as in Fig. 13. 

Fig. 16. — Model predictions for counts in three SIRTF bands (3.6//m, 24/xm and 70/xm). 
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Fig. 17.— Predictions of Model SI, Model S2 and Model S3 for color distributions of ISO 
sources. The detection limits in these bands are set to be /e./^m > 30 iiJy, /i5^m > 100 /iJy, 
/go/im > 100 mJy, fnoixm > 180 mJy, R < 24 mag, and K < 20 mag. The observed color 
distributions (histograms in corresponding plots) of ISOCAM 15/im sources in the HDF-N 
field are derived from data taken from Cohen et al. (2000) and Hogg et al. (2000). 

Fig. 18. — Predictions of Model SI, Model S2 and Model S3 for color distributions of 
SWIRE sources. The detection limits in these bands are set to be: froum > 6.15 mJy, 
f24nm > 0.45 mJy, fssum > 7.3 /iJy, fsum > 32.5 /iJy, R<25 mag, and K <20 mag. 

Fig. 19. — Panel a: 25/im local luminosity functions (LLPs) of three populations of dusty 
galaxies. Different hues denote different assumptions on the evolution of the sources: solid 
line — no evolution; dashed line — q=3.0; dashed-dot line — q=4.5. Panels b, c, and d: 
Plots of pairwise covariance of the LLP parameters. The eUipses show the 68% confidence 
intervals. The lines are the same as in panel a. 

Fig. 20. — The UV-B v.s. B-K correlation of galaxies. Data for ULIRGs (large crosses) are 
taken from Trentham et al. (1999). Other data are taken from the UV and IRAS selected 
sample of Xu and Buat (1996). 

Fig. 21. — The model SED of M82 compared to data. Data points: filled squares: collections 
of Silva et al. (1998); open squares: 2MASS and NED data. 
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